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Figure 1: Top: the experimental value o/sin^ Q \2 is compared with the predictions of exact Tri-Bimaximal 
(TB) or Golden Ratio (GR) or Bi-Maximal mixing (BM). The shift needed to bring the TB or the GR predic¬ 
tions to agree with the experimental value is small, numerically of order while it is larger, of order Ac 
for the BM case, where Xq = sin Oc with dc being the Cabibbo angle. Bottom: the experimental value of 
sin 013 in comparison with Xc or Xq. 


1. Introduction 

A long list of models have been formulated over the years to understand neutrino mixing. 
With time and the continuous improvement of the data most of the models have been discarded by 
experiment. But the surviving models still span a wide range going from a maximum of symmetry, 
as those with discrete non-abelian flavour groups (for reviews, see, for example. Refs. [ 1 , 2 , 3 ]), to 
the opposite extreme of Anarchy [ 4 , 5 , 6]. 

Among the models with a non trivial dynamical structure those based on discrete flavour 
groups were motivated by the fact that the data suggest some special mixing patterns as good 
first approximations like Tri-Bimaximal (TB) or Golden Ratio (GR) or Bi-Maximal (BM) mixing, 
for example. The corresponding mixing matrices all have sin^ 623 = 1 / 2 , sin^ @13 = 0 , values that 
are good approximations to the data (although less so since the most recent data), and differ by the 
value of the solar angle sin^ 612 (see Fig. 1 ). As the corresponding mixing matrices have the form 
of rotations with fixed special angles one is naturally led to discrete flavour groups. The relatively 
large measured value of sinOia has disfavoured TB and GR models because they in general predict 
sin 013 of the same order as the shift of the predicted sin^ 0i2 from the observed value, which shift 
is small for these mixing patterns. Instead in most models of BM the measured value of 0i3 ~ 9 ° 
[ 7 ] is natural. 

After the measurement of a relatively large value for 0i3 there has been an intense work to 
interpret the new data along different approaches and ideas. Examples are suitable modifications 
of the minimal TB models [8, 9 ], modified sequential dominance models [ 10 ], larger symmetries 
that already at LO lead to non vanishing 0i3 and non maximal 623 [ 11 ], smaller symmetries that 
leave more freedom [ 12 ], models where the flavour group and a generalised CP transformation are 
combined in a non trivial way [ 13 ] (other approaches to discrete symmetry and CP violation are 
found in Refs. [ 14 ]). 

Among discrete symmetry models, now that the value of sin^ 0i3 has been measured and found 
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to be not so small ([15]-[18]) the interest on BM mixing has been boosted [19, 20]. In this case 
the neutrino mixing matrix before diagonalization of charged leptons has sin^ d [2 and sin^ 623 both 
maximal and sin 0 i 3 = 0. The BM mixing matrix is in fact: 


Ubm — 


/-L__L 
^2 ^ 


2 

1 


2 

1 


0 ^ 
1 

^2 


V 2 2 V2 / 

The most general mass matrix corresponding to this mixing matrix is: 


mvBM = 



( 1 . 1 ) 


(1.2) 


where x,y and z are three complex numbers. One can think of a suitable symmetry or of some 
other dynamical ingredient that enforces BM mixing in the neutrino sector and that the necessary, 
rather large, corrective terms to 612 and 613 arise from the diagonalization of the charged lepton 
mass matrix. As well known, BM corrected by charged lepton diagonalization can be obtained 
in models based on the discrete symmetry S 4 [21, 22]. This idea is in line with the well-known 
empirical observation that 612 + 6 c ^ 7 r/ 4 , where 6 c is the Cabibbo angle, a relation known as 
quark-lepton complementarity [23]-[26]. Probably the exact complementarity relation becomes 
more plausible if replaced by 0 i 2 + &{ 6 c) ~ 7r/4 (which we call “weak” complementarity). In 
addition the measured value of 6\3 is itself of order 6 c'- 0 i 3 ~ 6 c j 3 / 2 . The relation of the neutrino 
mixing angles with 6 c could well arise in Grand Unified (GUT) models [26, 27] so that we will 
focus on GUT models in this work. 

In the present note we discuss two examples of GUT models of BM, one based on SU (5) and 
one on ^^(lO). The SU (5) model discussed here is more complete and indeed is based on a broken 
flavour symmetry that contains S 4 which imposes the BM structure in the neutrino sector and is then 
corrected by terms arising from the diagonalization of charged lepton masses. The present SU (5) 
model has the virtue that the quark mixing angles and the shifts from the BM values in the neutrino 
sector, all turn out to be naturally of the correct order of magnitude, expressed in terms of powers 
of Ac = sin 6 c> a’ la Wolfenstein, modulo coefficients of The ^^(lO) model is based on 

Type-II see-saw and the origin of BM before diagonalization of charged leptons is in this case left 
unspecified. Both GUT theories discussed in the following are variants of models already appeared 
in the literature, in particular by our group [22, 28]. We discuss the phenomenology of these 
models in the context of the present data and the comparison with other approaches like Anarchy 
and U{\)fn models. 

On the other hand, the relatively large measured value of 613 , close in size to the Cabibbo 
angle, and the indication that 623 is not maximal both go in the direction of models based on 
Anarchy, i.e. the idea that perhaps no symmetry is needed in the neutrino sector, only chance. The 
appeal of Anarchy is augmented if formulated in a SU{5)®U context with different Froggatt- 
Nielsen [29] charges only for the SU (5) tenplets (for example 10 ~ (a,Z 7 , 0 ), where a > 6 > 0 is 
the charge of the first generation, b of the second, zero of the third) while no charge differences 
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appear in the 5 (e. g. 5 ~ (0,0,0)). In fact, the observed fact that the up-quark mass hierarchies are 
more pronounced than for down-quarks and charged leptons is in agreement with this assignment. 
Indeed the embedding of Anarchy in the SU (5) (8> f/(l)fA? context allows to implement a parallel 
treatment of quarks and leptons. Note that implementing Anarchy and its variants in ^(^(lO) would 
be difficult due to the fact that all left-handed Standard Model (SM) fields of one generation belong 
to the same spinoral representation 16. In models with no see-saw, the 5 charges completely fix the 
hierarchies (or Anarchy, if the case) in the neutrino mass matrix. If RH neutrinos are added, they 
transform as SU (5) singlets and can in principle carry independent U{l)fN charges, which also, in 
the Anarchy case, must be all equal. With RH neutrinos the see-saw mechanism can take place and 
the resulting phenomenology is modified. 

The SU (5) generators act vertically inside one generation, whereas the U{l)fN charges differ 
horizontally from one generation to the other. If, for a given interaction vertex, the U{l)fN charges 
do not add to zero, the vertex is forbidden in the symmetric limit. However, the U{1)fn symmetry 
(that one can assume to be a gauge symmetry) is spontaneously broken by the VEVs v/ of a number 
of flavon fields with non-vanishing charge and GUT-scale masses. Then a forbidden coupling is 
rescued but is suppressed by powers of the small parameters X = Vf/M, with M a large mass, with 
the exponents larger for larger charge mismatch. Thus the charges fix the powers of X, hence the 
degree of suppression of all elements of mass matrices, while arbitrary coefficients kjj of order 1 in 
each entry of mass matrices are left unspecified (so that the number of order 1 parameters largely 
exceeds the number of observable quantities). A random selection of these kjj parameters leads to 
distributions of resulting values for the measurable quantities. For Anarchy the mass matrices in 
the neutrino sector (determined by the 5 and 1 charges) are totally random, while in the presence 
of unequal charges different entries carry different powers of the order parameter and thus suitable 
hierarchies are enforced for quarks and charged leptons by unequal tenplet charges. Within this 
framework there are many variants of models largely based on chance: fermion charges can all 
be non-negative with only negatively charged flavons, or there can be fermion charges of different 
signs with either flavons of both charges or only flavons of one charge. In Refs. [30, 31], given 
the new experimental results, a reappraisal of Anarchy and its variants within the SU (5) x f7(l)FN 
GUT framework was made. Based on the most recent data it is argued that the Anarchy ansatz is 
probably oversimplified and, in any case, nof compelling. In fact, suitable differences of U{l)fN 
charges, if also introduced within pentaplets and singlets, lead, with the same number of random 
parameters as for Anarchy, to distributions that are in better agreement with the data. 

2. A SUSY SU (5) model with S 4 discrete symmetry 

This model is a variant of those of Refs. [21, 22] to which we refer the reader for a detailed 
discussion and the technical details [32]-[35]. Here we concentrate on the differences, the moti¬ 
vations and the resulting phenomenology. This is a SUSY SU{5) model in 4-1-1 dimensions with 
a flavour symmetry S 4 ® ^ U{1)r iS> U (1)fa?, where U{1)r implements the R-symmetry while 

U(1)fa? is a Froggatt-Nielsen (FN) symmetry [29] that induces the hierarchies of fermion masses 
and mixings. The particle assignments are displayed in Tab. 1 . 

The formulation in 4-1-1 dimensions, with coordinate y in the fifth dimension compactified on a 
circle, allows a more efficient realization of GUT with less Higgs states, no doublet-triplet splitting 
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Table 1: Matter assignment of the model. The symbol br(bu) indicates that the corresponding fields live on 
the brane (bulk). 


problem and a better control of proton decay. In the present model it also introduces some extra 
hierarchy for some of the couplings. In fact, as indicated in the table, some of the particles are in the 
bulk (the first two generation tenplets T\ and T 2 and the Higgs and Hf) while all the other ones 
are on the brane at y = 0. Actually, to obtain the correct zero mode spectrum, one must introduce 
two copies, T\ 2 and T [2 with opposite orbifolding parity. The zero modes of Ti 2 are the quark 
doublets 2i 2 , while those of 7/2 are UI 2 and 2 - For economy of space only Ti 2 appear in table 
1. At leading order (LO) the S 4 symmetry is broken down to suitable different subgroups in the 
charged lepton sector and in the neutrino sector by the VEV’s of the flavons (pv, and The 

necessary proper alignment of the VEV’s is implemented, in a natural way, by the driving fields 
(Py, y/^, x^- The VEVs of the 6 and 6 ' fields break the EN symmetry. As a result, at EO the 
charged lepton masses are diagonal and exact BM is realized for neutrinos. Corrections to diagonal 
charged leptons and to exact BM are induced by vertices of higher dimension in the Eagrangian, 
suppressed by powers of a large scale A. While broken ^4 is at the origin of BM, U{\)fn, together 
with some higher dimensional effects, fixes the hierarchies of quark and charged lepton masses. 

With respect to Ref. [22] here we have modified fhe EN charges of fhe tenplefs from (3,1,0) 
fo (2,1,0). With this choice we aim at optimizing the neutrino mixing angles and to implement 
the "weak complementarity" relation at the expenses of a less accurate description of the first 
generation quark and charged lepton masses. We adopt the definitions: 


and 


A 



V(pv 


A ’ 

A 

A 

(e)_, 

(O') 

-t' 
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A 



s = 


s/tiRA 


<1 . 


where s is the volume suppression factor. 

It turns out that the simplest choice of setting all these parameters to be of (A h:)'- 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 


s = e = e' = t = t' = Xc, 


(2.4) 


leads to a good description of masses and mixings, as described below. Indeed by proceeding in 
exact analogy with Ref. [22], we have the following results. 
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2.1 Down quarks and charged lepton mass matrices 

For the down quark mass matrix, by only keeping the leading terms for each entry, one finds: 


md 


^aiiX^ ai2A^ £213^^ 




a33X 



(2.5) 


Here all matrix elements are multiplied by generic coefficients a,y of ^(1) with the exception of the 
(22) and (23) entries where the coefficients given by —c and c are equal and opposite. In the (32) 
entry the dots indicate that the lowest order potentially non vanishing is actually zero and the entry 
will get a contribution from still higher orders. For example, the entry (11) gets contributions from 
terms ~ X^). Thus the hierarchies arise from a combination of extra dimension factors, 

^4 and U{1)fn breaking. For the mass eigenvalues we have: 

nib-^v^^X'^, m, ( 2 . 6 ) 


For the charged lepton masses we have to take into account the introduction of the copies T/ 2 
of the first two tenplet fields, whose zero modes are differenl from fhose of ri ,2 and couple wifh fhe 
charged lepfons only. Therefore, all fhe operafors of fhe form FT^H^ have exacfly fhe same order 
1 coefficienfs whereas all ofhers confaining r /2 have differenl coefficienls. This Iranslafes in fhe 
following mass malrix: 

/a2iA^ a3iA^\ 

Me = ££ 12 ^^ —c'X^ . 1 Xv^ , ( 2 . 7 ) 

c'X^ £?33A j 

and fhe charged leplon masses are: 


, 0 t 2 


/ Wt = Wife ~ V^X 


run 


,0t 4 




nie 


,0i 6 




( 2 . 8 ) 


Nole lhal fhe b — T universalily is realized: here and in fhe following we obviously refer lo masses 
al fhe GUT scale as for example given in Ref. [36]. Nole lhal fhe predicled ralio mg/m^ ~ is nof 
perfecl, being foo large. 

The unilary lefl-handed rolalion Ud is oblained diagonalizing whereas fhe righl-handed 

one U£ is fhe charged leplon rolalion. Taking only fhe largesl conlribulion for each malrix elemenls, 
for Ui, which enters in fhe neufrino mixing malrix, Upnms = ulUy we have: 


^ 1 M 12 A MJ 3 A 

Ui ~ —MJ 2 A 1 0 

y-u\^X -ul2UljX^ 1 
so lhal 0 I 3 = 0 in Ibis approximation. 

2.2 Up quarks mass matrix 

Similarly Ihe symmelric up quark mass malrix is given by: 

^b\\X^ bi2X^ bi3X^ 
niu = bnX^ b22X'^ b23X'^ \ Xvl , 

yl7l3A^ b23X^ b33 


(2.9) 


( 2 . 10 ) 
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and the masses are: 

m, (2.11) 

The ratio ~ A v° /v° ~ A / tan j 8 implies that tan j 8 ~ 1 / A. In most of the cases the mass ratios 

are correctly reproduced, but, as already announced, there are some cases that need some moderate 
fine tuning. For example, we have ~ A^ which is too large, niu/nid ~ rndnis where, in 

reality the two sides differ by a factor of about 25. As we already mentioned we have chosen as 
a priority to fit the mixings rather than the masses. One can improve the agreement by somewhat 
relaxing the strict equalities in eq. (2.4). We now turn to describe the model predictions for the 
mixings first in the quark sector and then, after discussing the neutrino mass matrix, in the leptonic 
sector. 


2.3 The CKM matrix 


The CKM matrix is given by Vckm = U^Ud- The leading order expressions for and Ud 
obtained from the up and down quark mass matrices in eqs.(2.5) and (2.10) are of the form: 


/ 1 C12A M13A' 


and 




-c’|' 2 A 


1 0 

V-ct3A3 0 1 


( 2 . 12 ) 


Ud^ 


d\ 2 ^ \ 

1 d 2 -i^^ 

\{d\^dl^-d\^)X^ -d23X^ 1 ) 


-d*i2X 


(2.13) 


From these expressions we obtain the leading order form of the Vckm matrix with a pattern of the 
Wolfenstein type: 


/ 


Vckm • 


1 

—VJ 2 A 


V 12 A viaA^ 

1 V23A^ 


(2.14) 


\K2V23-'^13)^ -V 23 A 2 


1 


The Vij coefficients are related to the Cij and dij coefficients by: 

Vl2 = di2 —C*i 2 ', V 13 = (r/i3— 012 ^ 23 —C 13 ); V23=d23. 


(2.15) 


We see that for A = Ac the correct order of magnitude is derived for each Vckm matrix element 
modulo coefficients generically of order 1 . 


2.4 The neutrino masses and mixings 

The neutrino sector of the model is unchanged with respect to Ref. [22]. We therefore limit 
ourselves here to recall some important points. First note that in Tab.l there are no right-handed 
neutrinos. So the table refers to a model where the neutrino mass matrix is generated by the 
effective dimension 5 Weinberg operator. But a see-saw version is easily obtained by adding 3 right- 
handed neutrinos transforming under SU( 5 )x 5'4 x Z 3 as (l,3i, 1) and with charges f7(l)fi = -|-1 
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and U{1)fn = 0. The relevant phenomenology is quite similar, in particular the results on neutrino 
mixing. However, the neutrino mass spectrum turns out to be of moderate normal hierarchy type, 
with a LO lightest neutrino mass \mi\ larger than about 0.01 eV and, consequently, values of 
>3 — 4 10^^ eV. The deviation from a pure BM neutrino mass matrix is responsible for a 
softening of the lower bound on \mi \ and for an enlargement of the allowed values of the Ovj3j3 


rate. The experimental value of -y/r 




1/6 needs some fine tuning because in the model 


it should be of ^(1). This fine funing appears in mosf discrefe symmefry models because neufrinos 
musf be in friplef represenfafions of fhe discrefe group in order fo obfain BM or TB mixing efc. and 
fhen fhe mass eigenvalues are all of fhe same order of magnifude, barring cancellations. 

The neufrino mixing mafrix is obfained as Uy = uJUbm where Ue is given in eq. (2.9) and 
Ubm is the unitary matrix of BM. The results for the mixing angles are easily derived: 


sin 013 
sin^ 012 
sin^ 023 


\ui2 — Ui3\?i = |A| 

1 


1 

V 2 ' 

- - ^ Re{ui 2 + mi 3 )A = - -ReL 

1 

2 + ' 


The CP phase 5cp is 


5cp = 71 + arg{ui 2 - U 13 ) ^7i + arg{A). 


(2.16) 


(2.17) 


where the complex numbers A and £ are defined as A = — uij\X and 1,= ^ {ui 2 + Mia)^. 

The resulfs are graphically reported in Fig. 2 and compared wifh fhe experimenfal values of sin 0 i 3 
and sin^ 0 i 2 . We see fhaf, wifh X ~ Xc, the model realizes the "weak" complementarity relation and 
the experimental fact that sin 0 i 3 is of the same order than the shift of sin^ 0 i 2 from the BM value 
of 1/2, both of order Xc- It is interesting to observe that corrections to the BM pattern arising from 
next-to-leading order effects in the Yukawa couplings (higher order operators and shifts from LO 
flavon vevs) only affect sin^ 623 at the same ^{X^) as in eq.(2.16). We also see that, in general, the 
CP phase 5cp is not predicted, as the data only fix fhe absolufe value of A and nol ifs phase. If one 
could neglecf M 13 fhen A and £ would coincide. A very marginal agreemenf wifh fhe dafa would 
fhen demand fhaf bofh be aligned along fhe posifive real axis and, in fhis case 5cp = 7 i. 


3. Comparison with 5t/(5) 01/( 1) models 

If is inferesfing fo compare fhe previous model, which is rafher complicated involving SUSY 
SU (5), a non abelian S 4 symmefry and exfra dimensions, wifh a much simpler class of models based 
on SUSY SU(5) iS>U{1)fn [29]. As we have explicifly discussed a non see-saw version of fhe S 4 
model we will compare if wifh a non see-saw version of fhe U{1)fn models. In fhe following, only 
U{1)fn models wifh normal hierarchy are considered because, as shown in Ref. [37], U(l) models 
wifh inverse hierarchy (IH) fend fo favour a solar angle close fo maximal. 

In general we can label fhe U(l) charges as follows: 


10~ (tl,t2,0) 


5~(/i,/2,0) 


(3.1) 
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Figure 2: The complex numbers A = |mi 2 ~ M 13 |A and E = ^ (mi 2 + mi 3 )A determine sin 613 = |A| and 

sin^ 012 = 2 — ReT.. The CP phase is given by 5cp = 7Z + arg(A). The numbers shown approximately indicate 
the experimental central values. 


where 1, 2 refer to the first and second families. The Higgs field charges are taken as vanishing. 
There is also a flavon field 0 with charge -1 whose VEV (0) breaks U{1)fn- A set of charge values 
that lead to a good agreement with the observed masses and mixings are (this is the so-called 
case in [38]): 

10 ~ (3,2,0) 5 ~ (1,0,0). (3.2) 

The following mass matrices are obtained. For the up-type quarks: 

/ A^ A5 A3\ 

= A5 A4 A2 v« . (3.3) 

yA^ A^ 1 j 

Here A = (0)/A with A the large scale that suppresses the non rinormalizable interactions involving 
the field 0 and all entries are multiplied by coefficients that are complex numbers with absolute 
values of order 1. 

The down-type quarks and charged lepton mass matrices are one the transposed of the other 
and are given by: 

A^ A3 A3\ 

A3 A2 A2 v° . (3.4) 

Aliy 

Finally, the neutrino mass matrix transforming as 5 (8) 5 from the dimension 5 Weinberg oper¬ 
ator is given by: 
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A A 


My = 


A 1 1 
A 1 1 


A 


(3.5) 


We can now compare the two models starting from quarks and charged leptons. We observe 
that the structure of the mass matrices is very similar although not identical. In the S 4 model the 
suppression factors from the geometry in the extra dimension combine with those from the U{l)fN 
charges (which are different from those in the abelian model) to produce a similar pattern in the 
two cases. The mass matrices are not precisely the same in the two models but the predictions for 
the orders of magnitude of the mass ratios, expressed in terms of powers of A are identical. For 
both models we have in fact: 


mu 

md 

nis 

— 

rsj - 

rsj - 

Me 

Ms 

Mb 



me 



nif 


Mt: 

A^ 


(3.6) 

(3.7) 


Most of these orders of magnitude are correct if A ~ Ac but some, involving the first generation 
masses, are not, like niu/nic, mdlnis and nielmp,. Although the predictions for the mass ratios are 
the same, the ^4 model is superior, because its extra dimensional formulation solves the doublet 
triplet splitting problem and introduces corrections to the relation = mj of the Georgi-Jarlskog 
type. Also in the ^4 model m, / m,), ~ v:/(Av») tan j3 /A while the factor 1 /A is absent in the 
(7(1 )ca? model, so that only a moderate value of tanjS is needed in the S 4 model. The CKM quark 
mixing angles are also of the same order of magnitude in the two models and match the Wolfenstein 
pattern: 6^2 ~ A, ~ A^ and ~ A^. 

In conclusion, in the charged fermion sector the two models are rather comparable, with some 
advantages for the S 4 model. But where the latter is definitely superior is in the neutrino sector. As a 
result of the S 4 construction, the neutrino mixing pattern is dictated by BM corrected by terms from 
the diagonalization of charged leptons as detailed in eq.(2.16). The weak form of complementarity 
is realized as the shift of 6^2 is of the order of A ~ Ac and moreover also G^'j ~ Ac while @23 
deviates from the maximal value by terms of order A^. As already mentioned the observed value 

/ AtrP' 

of y/r = \ ~ 1/6 needs some fine tuning because in the model it should be of ^(1). In the 

y 

U{l)fN model the neutrino matrix is given in eq.(3.5). The diagonalization of charged leptons 
does not alter this pattern. For generic coefficients of ^(1) for each matrix entry, we would get 
that 0/3 ~ A, 023 ~ ^(i) which are good but also r ~ ^(1) and 0/2 ~ A which are bad. However, 
if, by accident, the 22 matrix element of my is of order A, then ^/r ~ A and 0/2 ~ i^(l): with a 
single fine tuning one fixes both problems. In any case, even accepting some amount of fine tuning, 
clearly there is no realization of weak complementarity. Finally in both models there are far more 
parameters than observables. This redundancy is less pronounced for the S 4 model but is still large. 


4. Bimaximal mixing in a 50(10) GUT model 

A challenging problem is that of formulating a natural model of Grand Unification based on 
50(10), leading not only to a good description of quark masses and mixing but also, in addition. 
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of charged lepton masses and neutrino mixing. In ^(^(lO) the main added difficulty with respect to 
SU (5) is clearly that all fermions in one generation belong to a single 16-dimensional representa¬ 
tion, so that one cannot separately play with the properties of the SU (5)-singlet right-handed neu¬ 
trinos in order to explain the striking difference between quark and neutrino mixing. A promising 
strategy in order to separate charged fermions and neutrinos in 50(10) is to assume the dominance 
of type-II see-saw [39] (with respect to type-I see-saw [40]) for the light neutrino mass matrix. 
If type-II seesaw is responsible for neutrino masses, then the neutrino mass matrix (proportional 
to) / (see eqs.(4. 1,4.3,4.4)) is separated from the dominant contributions to the charged fermion 
masses and can therefore show a completely different pattern. This is to be compared with the case 
of type-I see-saw where the neutrino mass matrix depends on the neutrino Dirac and Majorana 
matrices and, in 50(10), the relation with the charged fermion mass matrices is tighter. 

In Ref.[28], when the data suggested approximate TB mixing and a small value of sin^0i3, 
an 50(10) model has been studied based on type-II see-saw dominance. A detailed discussion 
of the general structure of this class of models can be found in the above article, together with 
a comparison with other approaches to 50(10) GUT’s. Here, given the relatively large value of 
sin 013 that has been recently measured, we reconsider this type of GUT model in the case of 
approximate BM mixing corrected by the charged lepton diagonalization. 

In renormalizable 50(10) models (a non necessary assumption only taken here for simplicity) 
the Higgs fields that contribute to fermion masses are in 10 (denoted by H), 126 (A) and 120 (Z). 
The Yukawa superpotential of this model is then given by: 

Wy = h^fx^H + f + (4.1) 

where the symbol y stands for the 16 dimensional representation of SO(IO) that includes all the 
fermion fields in one generation. The coupling matrices h and / are symmetric, while h' is anti¬ 
symmetric. The representations H and A have two SM doublets in each of them whereas £ has four 
such doublets. At the GUT scale Mqut, once the GUT and the B — L symmetry are broken, one 
linear combination of the up-type and one of down-type doublets remain almost massless whereas 
the remaining combinations acquire GUT scale masses. The electroweak symmetry is broken after 
the light Minimal Supersymmetric Standard Model (MSSM) doublets (to be called acquire 
vacuum expectation values (vevs) and they then generate the fermion masses. The resulting mass 
formulae for different fermion masses are given by (see, for example, [41]): 

Yu = h + r2f + r-ih\ (4.2) 

Yd = ri{h + f + h '), 

Ye = ri{h-3f + Ceh'), 

Yyo = h-3r2f + Cvh', 

where Ya are mass matrices divided by the electro-weak vev’s Vu^d and (a = 1,2,3 and c^, {b = 
e, v) are the mixing parameters which relate the Hu^d to the doublets in the various GUT multiplets. 

In generic 50(10) models of this type, the neutrino mass formula has a type-II and a type-I 
contribution: 

JYy = fvL-Mu^Ml, (4.3) 

fVR 
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where vl is the vev of the B — L = 2 triplet in the 126 Higgs field. Note that in general, the two 
contributions to neutrino mass depend on two different parameters, vl and vr, and it is possible 
to have a symmetry breaking pattern in ^(^(lO) such that the first contribution (the type-II term) 
dominates over the type-I term. The possible realisation of this dominance and its consistency with 
coupling unification has been studied in the literature [42, 43, 44, 45] and found tricky but not 
impossible [46]. The neutrino mass formula then becomes 

~ fvL ■ (4.4) 


Note that / is the same coupling matrix that appears in the charged fermion masses in eq. (4.2), up 
to factors from the Higgs mixings and the Clebsch-Gordan coefficients. Also note that the neutrino 
Dirac mass, proportional to Yyo in eq. (4.2), only enters in the neglected type-I see-saw terms and 
does not play a role in the following analysis. The equations (4.2) and (4.4) are the key relations in 
this approach. 

The 10 Yukawa couplings contributing to up, down and charged lepton masses in most models 
have a large 33 term, corresponding to the large third generation masses, while all other entries are 
smaller and lead by themselves to zero CKM mixing (because the 10 contributes equally to up and 
down mixing). Quark mixings arise from small corrections due to 126, the same Higgs represen¬ 
tation that determines / which in models with type-II see-saw is dominant in the neutrino sector, 
and to 120. Thus, in this approach, in the absence of 120, there is a strict relation between quark 
masses and mixings and the neutrino mass matrix. The presence of 120 dilutes this connection 
which however still remains important. In particular the deviations from BM mixing induced by 
the diagonalization of the charged lepton mass matrix, are typically of the same order as the largest 
quark mixing angle i.e the Cabibbo angle. An interesting question is to see to which extent the data 
are compatible with the constraints implied by this interconnected structure. 

For generic eigenvalues m,-, the most general matrix that is diagonalized by the BM unitary 
transformation is given by: 

f = UsMdiag{mi,m2,m3)ulM ■ (4-5) 

where Ubm is the BM mixing matrix given in eq.(l.l). In this convention Ubm is a real orthogonal 
matrix and all phases can be included in the eigenvalues m,-. Then the matrix / is symmetric with 
complex entries and, from eq. (4.5), one obtains (see eq.(1.2)): 


/ = 


(fi fi /i \ 

/i /o / 2-/0 ) 

\/i / 2-/0 /o / 


(4.6) 


with: mi=f 2 + V2fi , m 2 = /2 - \/2/i and m 3 = 2/o - / 2 . 

An important observation is that, for a generic neutrino mass matrix f, we can always go 
to a basis where f is diagonalized by the BM unitary transformation in eq. (1.1) and is of the 
form in eq. (4.6), in the same way as discussed in Ref. [28] for TB mixing. In fact, if we start 
from a complex symmetric matrix f not of that form, it is sufficient to diagonalize it by a unitary 
transformation U : = U^f'U and then take the matrix 

/ = Umf'diaA = U*BMU^f'UUl^. (4.7) 
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As a result the matrices / and f are related by a change of the charged lepton basis induced by the 
unitary matrix O = UU^j^ (in ^(^(lO) the matrix O rotates the whole fermion representations 16,). 
Since BM mixing is not a very good approximation to the data, in this basis substantial deviations 
from BM mixing must be generated by the diagonalization of charged leptons, with terms expected 
to be of ff{Xc)- At the same time also the quark mixings must be reproduced in agreement with the 
data. As the matrix elements of / enter both in the neutrino mass formula and in the corrections to 
the fermion mass matrices, this fact poses a non trivial problem of consistency, especially in view 
of the small values of the first generation masses. 

As one could decide to work in a basis where the matrix / is diagonalised by the TB matrix or 
by BM matrix (or in another suitable basis), this means that, for the measured set of data, the result 
of a best fit performed in one basis should lead to the same than the fit in other basis, because 
the only difference is that the set of parameters used in one fit are functions of the parameters of 
the other fit. So the x^ cannot decide whether TB or BM is a better starting point. However, since 
the first generation masses are very small some parameters must be precisely fine tuned in order to 
reproduce the small values of the masses. It is possible that one needs more fine tuning in one case 
than in the other. For a quantitative measure, in a given fit, of the amount of fine-tuning needed a 
parameter dpT was introduced in Ref. [28]. This adimensional quantity is obtained as the sum of 
the absolute values of the ratios between each parameter p, and its "error", defined, for this purpose, 
as the shift from the best fit value that changes the x'^ by one unit, with all other parameters fixed 
af their best fit values (this is not the en'or given by the fitting procedure because in that case all the 
parameters are varied at the same time and the correlations are taken into account): 

= (4.8) 

erri 

It is clear that dpr gives a rough idea of the amount of fine-tuning involved in the fit because if some 
\erri/pari\ are very small it means that it takes a minimal variation of the corresponding parameter 
to make a large difference on the x^- 

We report here on a comparative study of starting from / in the TB or in the BM basis. For 
the TB case the important difference with the detailed, complete discussion in Ref. [28] is that here 
we used updated experimental values for the neutrino mixing angles, in particular for sin^ 0 i 3 ~ 
0.022 ± 0.001, as most precisely measured by the Daya Bay experiment [17]. The result of a best 
fit performed in one basis should lead to the same x^ than the fit in another basis, because the only 
difference is that the set of parameters used in one fit are functions of the parameters of the other 
fit. So, as we have already stressed, the x^ cannot decide whether TB or BM is better. We have 
checked that the x^ is equal within uncertainties in the two cases, and this is true even for values 
of sin 013 somewhat different than the measured value, as can be seen in Fig. 3. However, since 
the first generation masses are very small some parameters must be precisely fine funed in order to 
reproduce the small values of the masses. It turns out that, for the physical value of sin^ 0i3, dpT is 
smaller in the TB case. A study of the fine tuning parameter when the fit is repeated with the same 
data except for sin^ 0 i 3 , which is moved from small to large, shows that the fine tuning increases 
(decreases) with sin0i3 for TB (BM), as shown in Fig. 4. 

A closer look at the details of the fine tuning parameter reveals that high dpp values are pre¬ 
dominantly driven by the smallness of the electron mass, combined with its extraordinary mea¬ 
surement precision. Moreover, due to the presence of mixing, the dpp coming from, for instance. 
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sin’013 


Figure 3: In the 50(10) model the resulting when starting in the TB or BM basis are equal within errors. 
Note that the minimum x^ value, X^ = 0.003, is obtained for sin^ 613 ^ 0.015, just a bit below the measured 
value sin^ 613 ~ 0.022. Nevertheless, as the minimum x} is quite shallow for sin^ 613 < 0.1, the fit does not 
exhibit any strongly preferred value 0 / 613 . 


the 33 component of h, which is mainly responsible for the top mass, is actually one of the largest 
contributions to the global dpT (due to its contribution to the electron mass) in both TB and BM 
scenarios. Although this might be surprising at a first glance, we emphasize that the dependence of 
the observables on the parameters is highly non trivial due to the off-diagonal elements of the mass 
matrices. 

In conclusion, as previously shown in Ref. [28], in this class of 50(10) models one can 
obtain a reasonable fit to the data. Then one can reinterpret the result as BM corrected by the 
charged lepton diagonalization and explicitly determine the corrective terms arising from the fit. 
However the model does not imply BM mixing as a starting approximation. In fact, one could as 
well focus on TB mixing and make a similar interpretation. To predict, before diagonalization of 
charged leptons, exact BM in the neutrino sector one would need additional dynamical ingredients. 
Independent of that, with the present value of sin^ 6 i 3 , a larger amount of fine tuning is needed 
in the BM case, as compared to the TB case, in order to reproduce the small values of the first 
generation masses. 

5. Summary and Conclusion 

Models of neutrino mixing based on discrete flavour groups have been extensively studied. 
After the recent measurement of sin^ 613 many of these models have been disfavoured, in partic¬ 
ular among those aiming at implementing TB mixing. But models based on 54 with BM mixing 
corrected by terms arising from the diagonalization of the charged lepton mass matrix remain as a 
viable and attractive possibility. In a GUT context, in these theories, it is also possible to implement 
the weak form of complementarity i.e. 612 -|- ^{Oc) ~ 7r/4 and to describe quark and lepton masses 
and mixings in a all comprehensive approach. Here we have discussed two examples of GUT mod- 
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Figure 4: In the 50(10) model the fine tuning parameter dpj increases (decreases) with sin^ 613 in the TB 
(BM) cases. For the physical value sin^ 613 ~ 0.022 it is about 4 times larger in the BM case. 


els of BM, one based on SU (5) and one on The SU (5) model discussed here indeed has 

a broken flavour symmetry that contains 54 and imposes the BM structure in the neutrino sector 
which is then corrected by terms arising from the diagonalization of charged lepton masses. The 
50(10) model is based on Type-II see-saw and the origin of BM before diagonalization of charged 
leptons is in this case left unspecified. We have discussed the phenomenology of these models in 
the context of the most recent data and their relative merits. We have then compared these models 
based on a large symmetries with models based on a minimum of symmetry where chance plays a 
central role, like Anarchy or models based on U{\)pi.j. The SU (5) model with broken 54 symmetry 
emerges as the most viable and predictive theory. 
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